Antishadowing and multiparticle production 



S. M. Troshin, N. E. Tyurin 

Institute for High Energy Physics, 
Protvino, Moscow Region, 142280, Russia 



Abstract 

We discuss role of absorbtion and antishadowing in particle production. 
Antishadowing, which leads to domination of elastic scattering at high ener- 
gies, appears to be consistent with growth of mean multiplicity in hadronic 
collisions. Moreover, we demonstrate possibility to reproduce power-like 
energy behavior of the mean multiplicity in the model with antishadowing 
and discuss physical implications of such behavior for the hadron structure. 
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Introduction 



Multiparticle production and global observables such as mean multiplicity and 
its energy dependence alongside with total, elastic and inelastic cross-sections 
provide us a clue to the mechanisms of confinement and hadronization. General 
principles are very important in the nonperturbative sector, in particular, unitarity 
which regulates the relative strength of elastic and inelastic processes. Unfortu- 
nately, there are no universal, generally accepted methods to implement unitarity 
in high energy scattering. Related problem of absorptive corrections and their sign 
has a long history (cf. [[1|]) and references therein). 

The choice of particular unitarization scheme is not completely a matter of 
taste. Long time ago the arguments based on analytical properties of the scattering 
amplitude were put forward [Q] in favor of the rational form of unitarization. It 
was shown that correct analytical properties of the scattering amplitude in the 
complex energy plane can be reproduced much easier in this form of unitarization 
compared to the most popular exponential form. Besides that rational form of 
unitarization leads naturally to prediction of the antishadow scattering mode [|3|]. 
Appearance of this mode is expected beyond the Tevatron maximum energy. 

Interest in unitarity and the corresponding limitations was stimulated under 
preparation of the experimental program at the LHC and the future plans to study 
soft interactions at the highest energies. Indeed, correct account for unitarity is 
also essential under theoretical estimates of the Higgs production cross-section 
via the diffractive mechanisms. The region of the LHC energies is the one where 
antishadow scattering mode is to be presented. It has been demonstrated that this 
mode can be revealed at the LHC directly measuring aei{s) and crtot{s) [0] and 
not only through the analysis of impact parameter distributions. Antishadowing 
leads to self-damping of the inelastic channels and dominating role of elastic 
scattering, i. e. crei{s) / atot{s) — ^ 1 at s — > cxd. Natural question arises about 
consistency of this mechanism with the growth with energy of mean multiplicity 
in hadronic collisions. Moreover, many models and experimental data suggest 
power dependence on energy of mean multiplicity^ and a priori it is not evident 
whether such dependence is compatible with antishadowing or not. 

In this note we apply the rational (fZ-matrix) unitarization approach ^ for 
consideration of the global features of multiparticle dynamics such as mean mul- 
tiplicity and role of absorptive correction. We show that it is possible to reproduce 
power-like energy behavior of the mean multiplicity in the model with antishad- 
owing and discuss its physical implications. 



Recent discussions of power-like energy dependence of the mean hadronic multiphcity and 
Hst of references to the older papers can be found in [|[ 
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1 Multiparticle production in the ?7-matrix approach 



The rational form of unitarization is based on the relativistic generalization of the 
Heitler equation of radiation dumping ^ . In this approach the elastic scattering 
amplitude satisfies unitarity equation since it is a solution of the following equa- 
tion 



F = U + iUDF (1) 

presented here in the operator form. Eq.[l] allows one to satisfy unitarity provided 
the inequality 

lmU{s, b)>0 (2) 

is fulfilled. The form of the amplitude in the impact parameter representation is 
the following: 

,x U{s,b) 
1 — iU[s, b) 

where U{s,b) is the generalized reaction matrix, which is considered as an in- 
put dynamical quantity similar to the eikonal function. Analogous form for the 
scattering amplitude was obtained by Feynman in his parton model of diffractive 
scattering [0]. 

In the impact parameter representation the unitarity equation rewritten for the 
elastic scattering amplitude f{s,b) at high energies has the form 

lmfis,b) = \fis,b)\' + r]is,b) (4) 

where the inelastic overlap function 

1 dcTinel 



Vis,b) = 



is the sum of all inelastic channel contributions. It can be expressed as a sum of 
n-particle production cross-sections at the given impact parameter 

Vi.s,b) = ^an{s,b), (5) 

n 

where 

1 d(y 

o-n{s,b) = ^-^5 ^n{s) = 871 bdban{s,b). 
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Inelastic overlap function is related to U{s,b) as follows 

lmU{s,b) 

Then the unitarity Eq. ^ points out that the elastic scattering amplitude at 
given impact parameter value is determined by the inelastic processes when the 
amplitude is a pure imaginary one. Eq. ^ imply the constraint |/(s, 6) | < 1 while 
the "black disk" limit presumes inequality |/(s, b)\ < 1/2 and the elastic ampli- 
tude satisfying this condition is a shadow of inelastic processes. The imaginary 
part of the generalized reaction matrix in its turn is the sum of inelastic channel 
contributions: 



ImU{s,b)=Y,Un{s,b), 



(7) 



where n runs over all inelastic states and 



Un{s,b) = / rfr„|f/„(s,6,{UP 



(8) 



and dVn is the n-particle element of the phase space volume. The functions 
Un{s, b, {^n}) are determined by dynamics of /ii + /i2 — > Xn processes, where 
{^n} Stands for the set of respective kinematical variables. Thus, the quantity 
lmU{s, b) itself is a shadow of the inelastic processes. However, unitarity leads to 
self-damping of the inelastic channels [|Sp and increase of the function \mU (s, b) 
results in decrease of the inelastic overlap function r]{s,b) when lmU{s, b) exceeds 
unity (cf. Fig. 1). 



1/4 



1/2 



lmf(s,b) 


unitarity limit 


\ black disk limit 


antishadow 


\ shadow 


scattering 


\scattering 



n (s,b) 


antishadow 


shadow \ 


scattering 


scattering \^ 



b=R(s) 



b=R(s) 



Figure 1 : Shadow and antishadow scattering regions 

Corresponding inclusive cross-section [T^ which takes into account uni- 
tarity in the direct channel has the form 

do „ r,.j,. '('•''•0 ,0, 
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The function I{s,b,^) is expressed via the functions Un{s,b,C., {^n-i}) deter- 
mined by the dynamics of the processes hi + h2 + 



n>3 

and 



})P (10) 



l{s,h,i)di = n{s,h)\mU{s,h). (11) 

The kinematical variables ^ (x and p^, for example) describe the state of the 
produced particle /i3 and the set of variables {^n-i} describe the system X„_i of 
n — 1 particles. 

Now we turn to the mean multiplicity and consider first the corresponding 
quantity in the impact parameter representation. The n-particle production cross- 
section cr„(s, 6) can be written as 

/ U\ Un{s,b) 



Then the probability 



anis,b) 



ainel{s,h) 



IS 



^"(^'^) = M7My ^^^^ 

Thus, we can observe the cancellation of unitarity corrections in the ratio of 
cross-sections cr„(s, h) and Oineiis-, h). Therefore the mean multiplicity in the im- 
pact parameter representation 



n 



[s,b) = ^nPn{s,b) 



does not affected by unitarity corrections and cannot therefore be proportional 
to?7(s,6). This conclusion is consistent with Eq. ([TTI). The above mentioned 
proportionality is a rather natural assumption in the framework of the geometrical 
models, but it is in conflict with the unitarization. Because of that the results of 
[ PTj ] based on such assumption and U -matrix unitarization should be taken with 
precautions. However, the above cancellation of unitarity corrections does not 
take place for the quantity n{s) which we address in the next section. 
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2 Growth of mean multiplicity 



As a starting point we use a quark model for the hadron scattering described in 
[p^. It is based on the ideas of chiral quark models. The picture of a hadron 
consisting of constituent quarks embedded into quark condensate implies that 
overlapping and interaction of peripheral clouds occur at the first stage of hadron 
interaction (Fig. 2). Nonlinear field couplings could transform then the kinetic 
energy to internal energy and mechanism of such transformations was discussed 
by Heisenberg [ [T3| ] and Carruthers Q]. As a result massive virtual quarks ap- 
pear in the overlapping region and some effective field is generated. Constituent 
quarks located in the central part of hadrons are supposed to scatter in a quasi- 
independent way by this effective field. 




Figure 2: Schematic view of initial stage of the hadron interaction. 

Massive virtual quarks play a role of scatterers for the valence quarks and 
their hadronization leads to the production of secondary particles. To estimate 
number of such quarks one could assume that part of hadron energy carried by 
the outer condensate clouds is being released in the overlap region to generate 
massive quarks. Then their number can be estimated by: 



N(s, b) oc ^ ^ ^""^ D^' ® (14) 

where mq - constituent quark mass, (kg) - average fraction of hadron energy car- 
ried by the constituent valence quarks. Function D'^ describes condensate distri- 
bution inside the hadron h, and b is an impact parameter of the colliding hadrons. 

Thus, N{s, b) quarks appear in addition to = n^^ + rih^ valence quarks. In 
elastic scattering those quarks are transient ones: they are transformed back into 
the condensates of the final hadrons. Calculation of elastic scattering amplitude 



has been performed in [ [T2| ] . 
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As it was already mentioned hadronization of massive N(s,b) quarks leads 
to formation of the multiparticle final states, i.e. production of the secondary 
particles. Remarkably, existence of the massive quark-antiquark matter in the 
stage preceding hadronization seems to be supported by the experimental data 
obtained at CERN SPS and RHIC (see [[T^] and references therein). 

Since the quarks are constituent, it is natural to expect a direct proportional- 
ity between the mean multiplicity of the secondary particles in impact parameter 
representation and number of constituent quarks appeared in the collision of the 
initial hadrons with given impact parameter: 

n{s,b) = aN{s,b), (15) 

with a constant factor a. The mean multiplicity n(s) can be calculated according 
to the formula 

^(^) = rri{s,b)r]{s,b)bdb ^ ^^^^ 
J^r]{s,b)bdb 

It is evident from Eq. (|16l) and Fig. 1 that the antishadow mode with the periph- 
eral profile of r]{s,b) suppress the region of small impact parameters and main 
contribution to the mean multiplicity is due to peripheral region of 6 ~ R(s). 

To make an explicit calculations we model for simplicity condensate distribu- 
tion by the exponential form, i.e. 



D,^~exp(-6/i?,). 
Then we have for the mean multiplicity 

n{s,b) = aii:iiMv^exp(-6/i?e). (17) 
The function U{s,b) is chosen as a product of the averaged quark amplitudes 

N 

U{s,b) = l[{fQ{s,b)) (18) 

Q=l 

in accordance with assumed quasi-independent nature of valence quark scattering. 
The 6-dependence of the function (fg) related to the quark formfactor FQ(q) has 
a simple form (/q) oc exp(— mgfe/^). Thus, the generalized reaction matrix (in a 
pure imaginary case) gets the following form [O] 



U{s,b) = ig 



1 + a- 



N 



exp(-M6/0, (19) 
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Figure 3: Energy dependence of mean multiplicity, theoretical curve is given by 
the equation n{s) = as^ (a = 2.328, 5 = 0.201); experimental data from [jT^]. 



where M = J2q=i ^Q- sufficiently high energies where increase of total cross- 
section is prominent we can neglect the energy independent term and rewrite the 
expression for U (s, b) as 

U{s, h) = ig {slmlY'^ expi-Mb/C). (20) 

After calculation of the integrals ( [T6| ) we arrive to the power-like dependence 
of the mean multiplicity n(s) at high energies 

n{s) ~ s^, (21) 

where 




We have two free parameters in the model, a and Rc, the freedom of their choice 
is translated to the free parameters a and 5. The value of parameter ,^ = 2 is 
fixed from the data on angular distributions JT^ ] and for the mass of constituent 
quark was taken the standard value mq = 0.35 GeV. From the comparison with 
experimental data (Fig. 3) on mean multiplicity we obtain that S has value 6 ~ 0.2, 
which corresponds to effective mass Mc = 1/Rc — O-Smq, i.e. Mc — m^r. 
It means that condensate distribution in the hadron is rather broad and does not 
coincide with the distribution of charged matter given by its formfactor. The value 
of mean multiplicity expected at the LHC maximum energy (^/s = 14 TeV) is 
about 110. Note that the numerical estimates for the total cross-section and the 
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ratio of elastic to total cross-section of pp-interaction at this energy in the model 
are the following: atot — 230 mb and (Jei{s) /atotis) — 0.67 0]. The latter value 
could help to detect antishadow scattering mode unambiguously. 



Conclusion 

It was shown that the model [ [T2[ ] based on accounting unitarity and extended to 
multiparticle production provides a reasonable description of the energy depen- 
dence of mean multiplicity leading to its power-like growth with a small exponent. 
This result is a combined effect of unitarity and existence of the phase preced- 
ing hadronization when massive quark-antiquark pairs are generated. It is worth 
noting again that power-like energy dependence of mean multiplicity appears in 
various models and is in good agreement with heavy-ion experimental data too 0. 

Multiplicity distribution P„(s, 6) and mean multiplicity n(s,b) in the impact 
parameter representation have no absorptive corrections, but antishadowing leads 
to suppression of particle production at small impact parameters and the main 
contribution to the integral mean multiplicity n(s) comes from the region of 6 ~ 
R{s). Of course, this prediction is valid for the energy range where antishadow 
scattering mode starts to develop (the quantitative analysis of the experimental 
data [jr|] gives the value: y/s^ ^ 2 TeV) and is therefore consistent with the 
"centrality" dependence of mean multiplicity observed at RHIC [[T9|]. 

In addition to the above conclusion, comparison with experimental data has 
shown that the peripheral condensate cloud of a hadron has rather large size. 

It is worth also noting that no limitations follow from the general principles 
of theory for the mean multiplicity, besides the well known one based on the 



energy conservation law. Having in mind relation ([T7|), we could say that the 
obtained power-like dependence which takes into account unitarity effects could 
be considered as a kind of a saturated upper bound for the mean multiplicity. 
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